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ABSTRACT 



We calculate the one-quark-loop amplitude for the low energy 77 vrvr collision in the 
context of the Nambu and Jona-Lasinio model with scalar and pseudoscalar four quark 
couplings to all orders in the external momenta. We show that the NJL predictions for the 
77 TT+TT" reaction are not far from the Born amplitude, which is close to the data, and is 
compatible with the chiral perturbation theory estimations. We determine the corrections 
given by the NJL model in leading order of 1/Nc to the chiral loop amplitude for 77 —>■ 
7r°7r°. Numerical results for the 77 tttt cross sections and for pion polarizabilities are 
given. 
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1 Introduction 



Currently there is much interest in calculations of the cross sections for the 77 — » tt+tt^ 
and 77 — >■ 7r°7r° modes. In these processes were considered in the framework of chiral 
perturbation theory (CHPT) 0-0] in the next-to-leading order. Moderate enhancement 
for the 7r+7r~ production near threshold and a small steadily rising cross section for the 
production were shown. For the charged channel the chiral calculation at the 



0^0 



TT TT 



next-to-leading order is in good agreement with experimental data For the neutral 
channel the theoretical picture is more complicated. Tree diagrams are absent at the 
level for 77 tt^tt^ which starts with meson-one-loop graphs. The one-loop result [Q] for 
^0^0 production deviates from the Crystal Ball data by 25-30% in the amplitude. Only 
the two-loop calculations (the order) ^ have reconciled CHPT with experiment and 
with dispersion calculations |@-||l3- Three new constants di which appear at this order 
m the L(6) CHPT Lagrangian have been estimated with resonance saturation. At present 
it is impossible to determine di from other processes. This still leaves an uncertainty in 
the CHPT calculation. 

Formerly an attempt was made to estimate the tree-level contribution which starts 



at order for Tr^Tr*^ production. Within specific models one can predict higher order 
constants of the effective meson Lagrangian at least at the tree level, i.e. to leading order 
in l/Nc- One can expect that tree-level contributions dominate over meson loops when 
they work at the same order in chiral counting. In this case, calculation of these tree-level 
contributions yields the main part of higher order corrections to the leading result 
for neutral pion production. Therefore there arises an opportunity not only to estimate 
the corrections appearing in the p^ approximation but also to advance in energy. This 
yields indirect information on the convergence rate of the chiral series. Besides, it is not 
impossible that in some cases strong compensations of contributions from box diagrams 
on the one hand and from a-exchange diagrams on the other can be observed in the p^ 
approximation. If so, the contributions of the order of p^ will be particularly important. 



In this paper, following in general the main idea of [|lT| we calculate the leading 1/iV, 



contribution to 77 Tin to all orders in external momenta and quark masses, in the 



context of the Nambu and Jona-Lasinio (N JL) model . A small number of parameters 
is an obvious advantage of the model as compared with the resonance saturation method 
where the number of parameters increases when new transitions and higher chiral orders 
are considered. Mathematical tools allowing the higher terms of the chiral expansion 



to be calculated are provided by the momentum-space bosonization method [13]. It is 



clear that other known methods ]|T5-]IT8| would lead us to the same results. We consider 



the NJL Lagrangian only with pseudoscalar and scalar couplings. The vector and axial- 
vector contributions are beyond the scope of this papeiQ. Also, the main bulk of events 
are peaked near on-shell values of the photons, so we present calculations only for real 
photons. 

In the recently published papers and ]|^ the NJL model is also used to investigate 



the process 77 7r°7r°. In contrast to this, we consider not only the neutral mode 
but also the charged one 77 vr+vr". Besides, in the above papers the amplitude was 



^We refer the readers to the papers and ||T^ where the TTTr-scattering amphtude has been calculated 
in the same framework using the NJL Lagrangian without spin-1 mesons and with them. The result clearly 
indicates only a small impact of a model extension like this on the final picture near threshold. 
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derived only to 0{p^) while we do not confine ourselves to this approximation. While 
analyzing the referee's comments on our work we learned of the results published in 
the preprint ||2^. Its authors employed another method to take into account the total 
momentum dependence of meson vertices |TB[ . They managed to solve the problem only 
numerically. In our paper we derive analytical expressions for the amplitudes in question, 
with full momentum dependence, and formulae for pion polarizabilities. This allows us 
to construct chiral expansions of the amplitude and to investigate the role of the first 
and approximations. 

The article is organized as follows. In Sec. 2 we set up the notation. In Sec. 3 we 
give a short description of the model and obtain general expressions for amplitudes of the 
77 — >■ TTTT reaction. In Sec. 4 we consider the specific neutral mode 77 vr^vr^, obtain the 
cross section for this process and calculate the polarizabilities of neutral pions. The same 
program is realized for the charged mode in Sec. 5. Finally, a summary and concluding 
remarks are presented in Sec. 6. 



2 Kinematics 

We consider the collision of two on-shell photons yielding two pions in the exit channel 

7(Pi, + 7(P2, e^.) 7r"(p3) + tt^Pa), 
where a, b are the isotopic indices. The matrix element for pion production is 

T{pi,P2,P3) = e^^^{PlY{P2)T^,u{Pl,P2,P?i), (1) 

where e is the electric charge. 

One can decompose T^^ into Lorentz and parity invariant amplitudes 

Tf,u = Cig^i, + C2P2fiPlu + C3P2f,P3u + C^PSfiPlu + C^P2,^,P3u (2) 

where the terms proportional to pi^ and p2i, are not considered since they drop out when 
contracted with the polarization vectors in (|l]), 

p''MPi)=P2^M=^- (3) 

From the Ward identities 

pIT^. = « = 0, (4) 

and the on-shell conditions pi = P2 = 0, P3 = pI = ^1 one gets the following constraints 
between the 

Cl + C2ipiP2) + CiipiPs) =0, Ci + C2ipiP2) + Csip2P3) = 0, 

Caipm) + c^ipm) = 0, ^4(^1^2) + C5(P2P3) = o. (5) 

Expressing all the scalar products (piPj) in terms of Mandelstam variables 

s={pi+P2Y, t=(pi-p3)2^ u={pi-p4y, s + t + u = 2ml, (6) 
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redefining C2 = A{s,t,u) and C5 = —sB{s,t,u), and with the following useful designa- 
tions, 

6« = ^ - nml, ^ = s,t,u, Y = ut-ml = uiti - sml, (7) 

one gets 

T^u = A{s, t, u)Ci^y + B{s, t, u)C2^,u, (8) 
with the gauge invariant tensors 

CT = p'2Pi - {g'^ (9) 

= - (^g'"" + + n.p^P'i + sp'.vi) . (10) 
The amplitudes {H^J) with equal (opposite) helicity photons are 

H^+ = -{A + mlB), H+^ = -B. (11) 

s 

The differential cross section for 77 tttt with averaged photon polarizations in the 
center-of-mass system are calculated by the standard expression 

/3(.) = (l-4m^/s)l/^ 
if(s,t) = + (12) 

The factor 5 is 1 (2) for tt+tt^ {tt^tt^). Let ^ be the scattering angle in the center-of-mass 
system. To compare our results with experiment, the maximum value of the integration 
variable (cos^)max is equal to 0.6 for ti^t^' [§] and 0.8 for 7r°7r° 0. 



3 The model 



3.1. The Lagrangian 

The original SU(2) flavor NJL Lagrangian [T^ with the fermionic degrees of freedom 
reinterpreted as the light up, u, and down, quarks incorporates the essential symmetries 
for the description of systems involving pions. We use it here in the form 

C = -M- + ^ [mf + (V5^75r^V')'] , (13) 

where G is the strong coupling constant of the four-fermion interaction and M. is the 
diagonal current quark mass matrix explicitly breaking the chiral symmetry of the La- 
grangian. We use = rhd = m for the matrix elements of M.. The electromagnetic 
interactions are introduced by the replacement — >• + ieQA^, where Q = {1 + 3r3)/6 
is the charge matrix in the isotopic space. 
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Pions appear as the Goldstone modes associated with the spontaneous breaking of 
chiral symmetry. The constituent quark masses (m) are generated by the gap equation 



m — TTi = 8mGIi. (14) 
We give all necessary definitions for quark-loop integrals Jj in the Appendix. We use the 



Pauli-Villars regularization introducing one Pauli-Villars regulator in such a way 

19 1 , that the basic scalar integrals Jj coincide with those of the usual regularization scheme 
with a covariant four-momentum cutoff, A. In the following we shall refer to the quark 
one-loops of N JL as tree contributions (in terms of meson fields) . To obtain the two-pion 
production amplitude on the basis of Lagrangian ( p!3[ ) one can use the momentum-space 
bosonization method or any other |]l4|-|jT8|. Let us stress that we calculate the tree 



level contributions at all orders in p^, in leading l/Nc- In chiral perturbation theory this 
corresponds to the coefficients of the operators in the tree level Lagrangian to all orders 
in a momentum expansion. 

The collective variables describing meson excitations are taken in the form of quark- 
antiquark field combinations linearly transformed with respect to the chiral group action. 
This method is completely equivalent to the existing nonlinear approaches but yields fewer 
Feynman diagrams in concrete calculations. We already discussed this issue in detail in 

Pf. 



3.2. The 77 — >■ tttt amplitude in the NJL model 

Our task is to calculate contributions from quark one-loop diagrams to the amplitudes 
A{s,t,u) and B{s,t,u) within the NJL model. Before going into detail, we make a few 
general comments. With the NJL model, one can calculate arbitrary A^-point functions 
to all orders in external momenta. These calculations correspond to allowing for all tree 
contributions, i.e. terms leading in l/Nc- The next step is to examine one-loop meson 



diagrams |2g] . However one has to confine oneself only to the first terms of the momentum 
expansion of the effective meson vertices derived from quark loops. Now it is not clear 
how to do it in a general case. Thus the following point of view seems reasonable. If 
there is a tree contribution, considering meson loop diagrams results in small corrections 
which can be ignored in the first approximation. And only if there is no tree contribution 
for whatever reason, consideration of meson loops is mandatory. This idea underlies for 



example the calculations in We shall follow it. 

The process 77 tt+tt" is known to have tree contributions in all orders of the 
chiral expansion. It is these contributions that we shall calculate by the NJL model 
here. For the neutral 77 — *■ tt'^tt^ channel, tree contributions appear only in the and 
higher approximations. At a level of the amplitudes A{s,t,u) and B{s,t,u) are only 
determined by meson loop diagrams calculated in fl], P]. This result does not depend on 
the model used and is of general character. We shall use it as a leading approximation]^. 

AL=-^'-^GAs), BZ=0. (15) 



•^We ignore the kaon loop contribution, since we work in SU(2); its contribution is small in the 
kincmatical region that we address M . 
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Here / is the value of pion decay constant /^r in the chiral hmit. The meson loop function 

Gtt (s) is given in the Appendix^. We use the notation of [0 for the loop functions. 
Further, when describing the process, we shall calculate quark loop corrections to the 
leading contribution. 

Let us turn to calculations. Contributing to the one-loop-order 77 mi amplitude in 
the NJL model are box diagrams of quarks, with two photons and two pions attached to 
the vertices, scalar exchange in the s-channel and pion exchange in the u and t channels. 
These diagrams are shown in fig.l. Pion exchange contributes only in the case of 77 

We proceed to calculate explicitly the contribution of the box diagrams, which is the 
most difficult to obtain. There are a total of six diagrams, with two distinct classes, one 
with photons being neighbors (four of this kind), see fig. la, and the one with photons 
alternating with pions (two of this kind). The complete amplitude for the box diagrams, 
T^t^^K is 

)}■ (16) 

Here g^, is the vrgg coupling constant. 



Al2{miy h{rnl) dp^ 

The factors /f* are the result of the trace over fiavor. They are equal to 



(17) 
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/f = y^afe, f^'=-{:^5az5,,-Uab). (18) 

As follows from these formulae, in the channels in question one has = = fi = 
10/9, f2 = —8/9, where we use the symbols N and C for the neutral and charged 
two-pion production respectively. 

The Lorentz tensors J^*j are the following two integrals, related to two topologically 
different box diagrams of fig. la 

'iVc\ f d^P 



J«(fci, h) = -I-jr[Siph,Sip + h^.Sip + h)^,S{p + ^3)75], (19) 



and 



jfX^i^ ^2, h) = Jy^ih - h, -k2, h - k2) = 

if) J A J^^^'^^^P^^^^^P + ^l)^55(p + k2h.S{p + ^3)75], (20) 

where the standard notation is used, 

Sip) = (21) 



*The index tt of the function Gtt (s) indicates that this function differs from the function G (s) given 
in the Appendix only by having the pion mass in place of the quark mass. 



6 



and iVc = 3 is the number of colors. As a consequence of Furry's theorem, the following 
identities hold 

JI^J{Pi,Pi+P2,P3) = 4m (P2,Pl +P2,P4) , 
JS){PuP1+P2,P4) = 4m (P2,Pl +P2,P3) , 
jjn){PuPl -P3,P4) = ji^J{P2,P2-p3,P4) ■ 

These identities were used to obtain (p!6D. 
One can find that 



j!n}iPuPl+P2,P3) + J!^J{Pi,Pi+P2,Pa) 



^Cif,uK2{s, t,u) + y {mlCi^y - €2^^^ ^3(5, t, 



u 



jfy{Pl^Pl -P3,Pa) = 



1 (s 
Here we use the following notation 



Y 



(22) 
(23) 
(24) 



(25) 



(26) 



Kl(s, t, u) 
K2{s,t,u) 

K3{s,t,u) 

^4(5, t, u) 



+ 



-g^^K,{s, t, u) + {p^,pl + ptpDi^th + Auh) + 

-{s- 2mlMp^ - p',p^,){Aj2 - AJ2), 
s 

khiO-himl)], 

hit) + hiu) + -{s- 2ml){uiAtl2 + hAj2), 
s 

3^3(5) + {S- 2ml)[Aj2 + AJ2 + l3{Pl,P3) + h{P2,P3) - 
Q4iPl,Pl+P2,P3) - Q4iPl,Pl+P2,P4)], 

shiPl, -P2) + tih{pi,p^) + uih{p2,P3) + {s- 2ml)Ke{s,t, 

- [th{pi,Pl+P2,P3) +Uh{pi,pi +P2,P4)] , 

(s - 2ml)[Atl2 + Auh + Qi{pi,Pi -P3,Pi) - h{P2,P3)] + 
2m^ r Y 

tihipi.p^) +Uih{p2,P3) - -^h{Pl,Pl -P3,P4) 



U 



s 

1 r 



hit) + h{u) - 2h{ml)\ + Qi{pi,pi -p3,Pi) - hiP2,P3), 

hit) + h{u) - 2/2 (m^)] + hiP3, -Pi) - Qi{pi,Pi +P2,P3) - 
1 



(27) 
(28) 
(29) 

(30) 
(31) 

(32) 
(33) 



Q4{P1,P1+P2,P4) 



2Y 



s - 2ml)[2tih{pi,p3) + 2uih{p2,P3) + 



+ shiPl, -P2) + shiP3, -Pi)] - 2ml[uihiPi,P3) + tlhiP2,P3) + 

+ sh{P3,-P4)]+ s[t'^h{PuPl+P2,P3) +u'^h{Pl,Pl+P2,PA)]} ■ (34) 

The rest of the notation is given in the Appendix. Note again that we systematically use 
the conditions (0) when deriving ( p5[ ) and 
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Now let us evaluate the scalar exchange, which contributes only to the s-channel, see 
fig. lb. Here one can obtain 



rp((7) 
fJ.U 



[2/2(5) + is- 2ml)h{ps, -p,)] 



+ {Am? 



(35) 



The function Q-i{s) results from calculation of a triangular quark diagram describing two- 
photon decay of a scalar meson. It is given in the Appendix. The expression in the 
numerator is well known. It corresponds to the crTrvr vertex. The denominator stems from 
the scalar particle propagator. It is equal to [rn?^{s) — s]/4(7^(s). The derived expression 
(35) is gauge-invariant. 

Finally, we evaluate pion exchange, , which contributes to the t and u channels in 
the case of production of charged pions, see fig.lc. For on-shell photons we obtain 



rp(TT) 



+ 



Qj2nr—l/fah £ah\ 

°9-K^-K l/i - J2 )\ 72 



th{t) ~ mlh{ml) 



+ 



ul2{u) -mj2{'mj 



(36) 



In eq. (^) the pion propagator is cancelled by one of the ■jTiTr vertices. 

The sum of (|16D and (|36D meets the general requirement of gauge invariance because 
it is proportional only to the Lorentz tensors Ci'^ and 



with 



1 

Y' 



-2/€7£r 



s 



Y 



Apart from the known functions K,i{s,t,u) (see (29)- 



K7 
K4 



s 



hit) , hiu) 



Ml 



1 + 



£^'')(s - 2m^)K5. 
, we also use new ones 

2Mltl) 



rntis 



StiUi 



K2 - siAth + Kh] 
- siAth + AJ2) 



(37) 

(38) 
(39) 

(40) 

(41) 
(42) 



Thus, we have all that is necessary for exploring physical consequences of the process 
under consideration. Before doing this, we point out that besides direct computation we 
also use another method for calculating amplitudes, namely projection of the sought-for 
amplitude T^^{s,t,u) by means of convolutions with tensors Ci'^ and £2'^. A symbohc 
algebra program has been written for this purpose in REDUCE. First, the contributions 
to the amplitudes are calculated as functions of tensor integrals for all diagrams. Then the 
well known reduction method |2^ is applied in order to transform the tensor integrals to a 
combination of a few basic scalar integrals (^)-(|y). Finally, the amplitudes thus obtained 
are projected onto the gauge invariant tensors (j|) and (|IOD. The last step is necessary in 
order to have few compact and simple expressions, which are independent. This represents 
a consistency check for the methods used and in particular for the regularization scheme 
applied to the quark loop integrals. 
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4 The 77 n n mode 



4.1. Chiral expansion of the amplitude 



Theoretically of greatest interest is the neutral mode. Here the leading approxima- 
tion is described by meson loops (see (0)). Yet, as was shown in the literature, this is 
not enough to describe the process quantitatively even in the near-threshold region. Let 
us elucidate the role of tree corrections. To this end, we isolate terms of order p^ and p^ 
from our general result. These terms are most important near the threshold. We start by 
expanding the amplitudes in a series of momenta. For box-type diagrams we get 

"^^^^ 4^-(™-|; 



rp{box)N 



+ 



h+ 



+ 



4m: 



1. 

2 



' s — lAml 



+ ... (43) 



where the superscript stands for neutral channel. From here on we use the following 
notation 



167r2 



3hi 



h = HO) = ^ ( In 
h = h{0) 

^4 = ^4(0) = -^, h- 

^ ~ 3207r2m6(A2 + m2)2 
From the scalar exchange we have 



A2 



A2 + m^y 
,2 \2 



A^ 



^A2 + 3m2' 
A2 + m? 



hi, 



(44) 

(45) 
(46) 
(47) 



rp{o-)N 



C 



2\2 



+ 



9Z, 

(s — m,^^ 

24m4 
7s{s — ml 

240m2 



2 7s 
3^^ " 60^^ + 



s — ml ^ s — 2mi 



6m2 



-/3 + 



/3 + 



3s2 — llsmi + 8m^ 
36m2/2 



Il + 



3/2 

sis 



2ml 



9/1 



-'3 



'2 

40m^ + 14sm2-17s2 
"-'4 H 77:777 hh 77rh + • • • 



120/9 



42 



(4J 



As might be expected, by summing (|43| ) and ( ]48| ) the constant terms (2/3/3) are cancelled 
(absence of tree contributions at a level of p^), and others are grouped into a combination 



rp(^box+a)N 

fJ,U 



40^ 

9Z^ 



C 



(s-mt 



1 + 



s — mz 



+ 



+ 



2mi 



Qm^ \ 4m2 
3s2 — 1 Isml + 8m^ \ 



^ 2 7s(s — m2)\ 
■4m2 + \^ /4+ 



60m2 



12m2 

40ml + 14sm2 -17^2 



7^2 

3/2 



+ 



s[s- 



9/2 



2m^) 



1. 

2 



120/2 



/.- 



/3/4 + ^ (.^ 



3Miti + 8ml /; 



14m2 



+ . . . 



(49) 
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Now let us switch to the chiral expansion of the obtained amphtude. This is an expansion 
in external momenta and masses of current quarks. We discussed the problem in detail 



within the NJL model in |T^, where all necessary formulae can be found. The quantities 
corresponding to the leading terms of the chiral expansion will be labelled by an "o" above 
the respective symbol (except the chiral limit for the constant /tt traditionally designated 
by /). By algebraic transformations, we find (see (H)) for the amplitudes A{s,t,u) and 
B{s, t, u) 



A^{s,t,u) 



o 2 

40 m 



+ 



+ 



+ 




1 - 



om 



+ 



o 2 o 2 

3m 



4m 



o 2 o 2 o 2 

7s{s-m^) 2m^(s-3m^) 



(A2+m^)(A2 + 3m\ 



+ 



o 2 

60m 



o 2 

m 



+ 



o3 
o2 



o 2 



o 4^ 



6sm„ +12m. 



o o 

hh 



o 2 



o 2, 



120 h 

o 2o 

20 m 1 4 
9P~ 



16m^ (9s - 20m^) - 17s 



o 2 o2o2o2 o2o o2 

m^ 3m mh^ m^/^ 14m^— s o 

TV TT'^l TT-' -J I TV Y 

o 2 o o ~r O 

m A4/i2 3/2 3/4 



/5 



(50) 
(51) 



These expressions describe the tree-level contributions to the amplitudes A'^{s,t,u) and 
i?^(s, t, u), which involve terms of orders and p^. At A ^ 00 the term proportional to 
1 in the amplitude -B^(s, t, u) coincides with the expression derived in |Tl|. The terms of 
order p^ in the amplitude v4^(s, t, u) differ from the result of [11]. The discrepancy is due 
to the contribution from the a meson exchange diagram, which was not considered there. 
In the p^ approximation the CHPT Lagrangian has three constants di, 62 and ^3 



entering in the description of the process 77 
order p^ in the process in question has the form 



0. In terms of di the amplitude of 



A(6)(s,t,M) 



40 

80rfi 
W 



o 2, 



[s - 2m^){Ad2 -di- 84) + 2{s-m^){di + Ad 



(52) 
(53) 



It is easy to see that the NJL model yields the following values for the constants 



di 



do 



dp. 



2fh2 



[IQnm)'^ 
P 



-1.26 X 10 



-4 



(167rm) 



hi +hi 



o 2o2' 

3m hi 



Phi 



2(167rm)2 



2.89 X 10" 



0.96 X 10" 



(54) 
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The analytical expressions for these coefficients fully agree with the corresponding ex- 
pressions derived in |2^, if one sets qa = ^ there. For quantitative evaluation we used 
the numerical values of the model parameters fixed earlier in |T3|: m= 221.2 MeV, / = 
88.6 MeV, 'nn^= 141.5 MeV and A = 1 GeV. Note that these values correspond to 
the leading terms of the chiral expansion of the physical quantities chosen as follows: 
= 139 MeV, U = 93.1 MeV, m = 241.8 MeV, A = 1 GeV. 
An equivalent (ij-related set of parameters is also used. It includes constants ai, 02 
and h. 



ai 



0-2 

;i67r2/2)2 

h 



(167r2/2)2 
The values we obtain for them are 



a2 



10(47r/)- 

o 2 

9m 
10(271/)- 

o 2 

9m 



5(27r/)- 

o 2 

9m 




o 2o2 

3m hi 
27r2/2 



hi ~2h2 



4/ii +3h2 



o 2o2' 

3m hi 



(55) 



-59.5 



35.6 



h2= 3.5 



(56) 



Let us compare these values with the results of where the parameters ai,a2 and h 
were calculated within the extended NJL model. 



(ENJL) 

(ENJL) 
"2 

^{ENJL) 



-23.3, 

14.0, 

1.66. 



(57) 



It is quite simple to account for the discrepancy. In the extended NJL model Oi, 02 contain 
a factor g\ ~ 0.37 and h contains a factor qa — 0.61. This factor arises from allowing 
for vector mesons. In the NJL without vector mesons it is equal to 1. This is the reason 
for larger values of the parameters in (^6]). Yet, this does not mean that the NJL model 
without vector mesons cannot be used to describe the process in question qualitatively. 
This only indicates that chiral expansions in the NJL model converge slower than in the 



ENJL model, which we already know from [IS]. Really, this means that if we want to 
compare the NJL result with the CHPT calculations, where chiral expansions converge 
quite fast, we must compare the total NJL result, which takes into account all orders of 
the chiral expansion. We shall return to this matter again when discussing numerical 
values of pion polarizabilities. 



4.2. Polarizability of neutral pions 
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To obtain the pion polarizabilities we go to the crossed channel, Compton photon-pion 
scattering. This is easily done by trivial replacement of variables. The general expressions 
for the polarizabilities are obtained as 



PI 



2m, 



s,4i,ui=0 



(58) 

[oL^ + = -am^5*|^_t^,„j=o (59) 

with [i = N,C) and the coupling constant a = e^/(47r) = 1/137. In the case of charged 
pions, (C) calculated below in 5.2, the Born term is removed from the amplitudes. Dif- 
ferent phase conventions are used in the literature. We adopt the one chosen in This 
is the Condon-Shortley phase convention. The one-loop result is 

a 



0, 



4.0. 



From here on we express the numerical values of polarizabilities in units 10~^ fm'^ 
this result one must add the and p^ tree contributions obtained by us. 

o 2 



(60) 
(61) 

To 



N 



10am 



+ 



N 



+ 



+ 




o 2 o 2 

3 m m^ 



(A2+m^)(A2 + 3m^) 



+ 



1.14 



(62) 



o 2o2 

3m hi 



5 m, 



Sm 



O 2 o 2 

3 m m^ 



2 o , . , „ o 2 o 



o 2 o 2o3 

15 m^m hi 
(2vr/)4 



0.92 



+ 



(63) 



15m (A2+m )2 

The numerical evaluation results are given in Table 1. It is evident from the table how 
important it is to calculate polarizabilities of neutral pions with allowance made for all 
terms of the chiral series and not only for its first terms. This can be most clearly 
demonstrated by examining the sum and the difference of polarizabilities. The pion loop 
is known not to contribute to the sum (ott + Ptt). In other words, this quantity stems only 
from many-loop contributions (meson two-loop, three-loop contributions, etc.,). 



(a^r + Pn)cHPT ~ X! 



n=2 



(4^ 



E 



(64) 



where the index n denotes a contribution from n-loop meson diagrams. The quantity 



is a sum of the contact contribution s*" 
and the proper loop contribution s^^^ 



(n) 
(loops) _ 



'(2) 



= h[_ in terms of notation from [0) 



5(n) 



r I (loops) 
S(n) + S(n) ■ 



(65) 
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Table 1: Neutral pion polarizabilities in the NJL model in units of 10 ^ fm'^. The results are 
compared with CHPT. 





1 loop 


tree 


total 


CHPT 










all others 




[7] 







1.44 


-0.30 


0.05 


1.19 


1.17 ±0.30 




-1.0 


-1.63 


0.92 


-0.45 


-2.16 


-1.90 ±0.20 




-0.5 


-0.09 


0.31 


-0.20 


-0.48 


-0.35 ±0.10 




0.5 


1.54 


-0.62 


0.25 


1.67 


1.50 ±0.20 



In they calculated the first term of (|6^) S{2) = 1.0 ±0.17. As in (|65|), here the first term 
corresponds to the contact contribution and the second term to the loop contribution. The 
total tree contribution, which is calculated in the NJL model, corresponds to the sum of 
the contact terms in a similar CHPT calculation. For the sum of polarizabilities of neutral 
pions we obtain the value 



oo 



K + /5.)Ll = 1-19 ~ E^(n)- (66) 

n=2 

Hence it can be concluded that the resonance saturation hypothesis, which is used to 
determine the contact contributions in CHPT and results, among other things in S^2) — 
1.0, holds good here because it leads, as it must, to fast convergence of the CHPT chiral 
series. It follows from our analysis that the sum of all other terms of this series accounts 
for only 20% of the leading contribution. Unlike the CHPT, where already the first term 
of the chiral series practically determines the entire result, the NJL model without vector 
mesons allows this at the second step. Moreover, if we confine ourselves to the first two 
terms of the chiral series in calculation of pion polarizabilities by the NJL model, we shall 
fully reproduce the picture we have in the CHPT with and fixed in conformity with 
the resonance saturation hypothesis. The parameter is not associated with the pion 
polarizability, so we cannot get any more severe constraints on its value as compared with 
the already known ENJL result. 

An interesting picture emerges for the difference of polarizabilities (ott — Pn)- Here we 
have 

{an - Pn)cHPT = ^(n) = 77^^ d{n)- (67) 
n=l l^^JJ „=1 

As in (|64D, the quantity d{n) is a sum of the contact contribution d^^^ [(T^^) = ±86^ = h\_ 
in terms of notation from [0]) and the proper loop contribution In they obtained 

D^2) = —0.58 and D^^"^^^ = —0.31. In the NJL model, if we again confine ourselves to the 
sum of the first two terms of the chiral series, we get a close value [{D^^)) ± {D^3))]njl = 
—0.70 ~ {D^^2))cHPT- However it is evident from Table 1 that the chiral series of the NJL 
model for the difference of polarizabilities (a^ — (3-,^) converges noticeably worse. In this 

case we get 

oo 

(«.-/5.)j^jL = -l-16 ~ E^(n), (68) 

n=2 

which is twice as large (in absolute value) as the leading (-D[2)) result of the CHPT. Hence 
it follows that in this case the resonance saturation hypothesis is less successful. This is 
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no surprise. The dispersion calculation methods leading to the sum rule for (a^ + /?^) 
are known to yield quite an exact result for this quantity. This is not the case with 
the difference (a^r — /^vr) which features uncertainties in determination of the asymptotic 
contribution [^. It is in these cases that the complete NJL model calculations are of 
particular interest. As we have just shown, they allow quite concrete conclusions to be 
drawn about the efficiency of the approximations used. 

Let us also stress that the correction to the electric polarizability at order is very 
small owing to strong compensation of the terms related to a exchange (see the expression 
in parentheses) 



2o 



A(^%f = ^ - 1 = -0.09. (69) 

(127r/m)2 \ 2ti^P ) 

Therefore and higher order terms turn out to be essential here. Our calculations show 
that the contribution of these terms compensates almost fully the result. It indicates 
that the two-loop meson graphs are important for this value. 

The process 77 — >■ tt^tt^ has been used to obtain information on the polarizabilities. 
We give here the results obtained from the Crystal Ball data by Kaloshin, Persikov and 
Serebryakov. 

{a^ + (3^)^ = 1.14 ± 0.08 ± 0.16, |2§ 
{a^-^^f = -1.1 ±1.7. H 

4-3. Interaction cross section 

Before showing the results, we need to explain the parameters used. The values of 
the parameters m, change at each order of the chiral expansion, but can always 



be expressed in terms of the leading chiral ones, m, /, [|^, which has been done in 



expressions (|^)-(|^ and in equations (|7T])-(^) below. The Mandelstam variables are 
always constrained by s + ti + Mi = and since the pion mass changes order by order, so 
does also the kinematically allowed region for s, t and u. For that reason the pion mass 
which appears in the polarization tensors C'l'^ and £2'^ is the one calculated at a specific 
order. The same is true also for the phase space factor ([l^). In this way we get the 
same result by summing an infinite series of the chiral expansion or calculating the full 
amplitude. The pion loop has been calculated only at leading order and we have to take 
the pion mass of the specific tree order to which it is added. In Fig. 2 we display the data 
for the cross section a{s, \ cos6\ < 0.8) as determined by the Crystal Ball collaboration 0. 
They are shown as a function of the center-of-mass energy W = ^/s. The dotted line is a 



well-known result arising from consideration of the leading pion loop (15). The dot-and- 
dash line is the result at which we arrive on taking into account p^ terms. The dashed 
line shows the data derived from the previous curve by including order p^ corrections 
(|5l|)- (|5l|) . The solid line shows the full result. The tree contributions are clearly seen 
to be responsible for enhancement of the scattering cross section in the near-threshold 
region, the tree p^ approximation accounting for almost the entire effect. Diagrams with 
0" exchange dominate among the diagrams contributing to this result. Consideration of 
higher-order tree corrections practically do not affect the cross section in the low-energy 
region. The different thresholds correspond to the different values of the pion mass at a 
certain order of the chiral expansion. 
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5 The 77 tt^tt mode 



5.1. The result up to 

Let us begin with the leading-order result. Here the NJL model yields 

hi . 4 



A 



c 



+ 



B 



+ 



{2Tcfy uiti 

The ellipses indicate higher contributions. The amplitude resulting from insertion of 
B'^ in (H) is the Born term and of order p^. The term A'-^ leads to the p'^ correction 
of the amplitude. Our task is to calculate corrections to these expressions arising from 
consideration of NJL tree diagrams. They can be derived from the general result given 
in Sec. 3. Here we shall dwell upon the role of the p^ and p^ corrections which account 
for the major part of the additional contribution. 



5.2. The result up to p^ and higher orders in the NJL model 



Let us examine corrections to the leading terms of the chiral expansion of the 77 
TT+TT" amplitude. Here we have 

o 2 



A^{s,t,u) 



+ 



o 2 

8m 

w 

O 

/3 



/3 1 + 



5 s- 



o 2, 



18m 



_ (5.-13^:)+^ 

9/2 6/3 



B^isA.u) 









/ 2 \ 2- 




1 H 


\ 


3m^/ii \ 


1 ^^1^1 










10(67r/m) 



(71) 



(72) 



but 



One obtains the polarizabilities for the charged pions from equations (^) and 
with the Born term removed from the amplitude. We derive the p'^ and p^ orders to the 
polarizabilities of charged pions to be 

c _ 2a/ii I . 5m! 11 ^ / 



m^ (47r/)2 



1 



o 2 

18m 



(a^ + 



11 
20 



— + 2 



o 2n 

3m 



o 2o 

3m hi 



o 2 

A2+m 



27r2/2 



-4.57, 



7a m„h2 



0.41. 



(73) 



(74) 



10(67r/m)2 

Generally, it is these contributions that determine the amount of the electric and magnetic 
polarizabilities of charged pions. Table 2 presents numerical evaluation results for the 
polarizabilities. We also included corrections from p^ terms in the table. Analytical 
expressions are quite cumbersome, so we do not show the whole of them. An exception 
is the sum (a^ + (3.,^)'" , for which the order p^ contribution is a leading one and the first 
correction at order p^ amounts to 25%. The NJL model yields for it 

2 



7a m^h2 
10(67r/m)' 



125 m^ 

1 J + 

168m 



o 2 o 

(2vr/)2 



o 2 o 2 o2 

13m m^hi 
7A4 h2 



0.31. 



(75) 
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Table 2: The charged pion polarizabihties in the NJL model in units of 10 ^ fm^. The results 
are compared with phenomenological data. 





tree 


total 


exp. 








all others 







0.41 


-0.10 


0.02 


0.33 


1.4 ±3.1 ±2.5 [31 
0.33 ±0.06 ±0.01 1 
0.22 ±0.07 ±0.04 1 


29] 
291 




11.62 


-2.07 


0.82 


-0.36 


10.0 


15.6 ±6.4 ±4.4 [P 
4.8 ± 1.0 |10| 




5.81 


-0.83 


0.36 


-0.17 


5.17 






-5.81 


1.24 


-0.46 


0.19 


-4.84 





3TI1 from the reaction tt A 



Apart from the known results obtained in Serpukhov 
'n'~'yA, the experimental body of the table incorporates the results of analyzing angular 
distributions of pions in the reaction 77 — vr+vr" (CELLO and MARK-II data) within 
the framework of the unitary model for helicity 2 amplitude ^9]. As seen from Table 
2, the results of numerical calculations agree with the available experimental data. The 
electric and magnetic polarizabihties are mainly determined already at order p^. Only 
the sum of these quantities, which equals zero at this order, is determined by higher 
and corrections. 



5.3. The amplitudes: numerical results 



Let us turn to the diagrams illustrating the general formulae derived earlier. Figure 3 
displays the charged pion production cross section as a function of energy in the center- 
of-mass system of interacting particles. Integration over the scattering angle covers the 
interval | cos^| < 0.6. We show four curves. The dotted curve corresponds to the Born 
term derived from in formula ([701). Taking into account corrections up to p^, ([70|) 
and up to p^ (see formulae in (|7T1) -([72|)) we get the dot-and-dash and the dashed curves 
respectively. Full calculations are represented by the solid curve. As was expected, the 
cross section for production of charged pion pairs features some enhancement caused by 
p^ terms. The total result is close to p^ estimations, in the low energy range. For higher 
values of the energy the scalar resonance starts to show up in the higher orders in momenta 
and in the full result 



6 Conclusions 



To summarize, we have calculated the quark one-loop contributions to the 77 —>■ iriv 
processes, using the SU(2) flavor NJL model, with scalar-isoscalar and pseudoscalar- 

^This behavior differs from the result obtained in where the scalar resonance was hidden. This 
difference is a consequence of the different regularization schemes used. In order to have the right low 
energy behavior in the Compton photon-pion processes considered here one has to use Pauli-Villars 
regularization. 
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isovector four-quark interactions. These calculations correspond to full consideration of 
tree contributions in the examination of the above processes. We present the amplitudes 
in a compact and analytical form which is used to compute the associated cross sections 
and the polarizabilities of pions. 

The chiral expansion of the 77 — > vr^vr^ amplitude starts with terms. There are 
no tree contributions in this approximation, therefore the leading contribution is that 
from one-loop pion diagrams. We calculated the highest tree corrections appearing at 
order and showed that in the near-threshold region the cross section for this process 
can be satisfactorily described by one-loop pion diagrams (the order) with the tree p^ 
contribution. The conclusion as to the smallness of the highest contributions (beginning 
with terms) to the cross section of the process under study in its near-threshold region 
is of interest to us. It is this kind of evaluation that the quark models are well suited 
for. For example, careful examination of this problem within the CHPT at order p^ alone 
would require complicated calculations of the three-loop meson diagrams. 

The NJL model predicts some enhancement of the cross section for 77 tt+tt^ near 
threshold. This fact is in good agreement with the known CHPT result [0]. The increase 
of the cross sections in the full result for higher energies is an artifact of the model. 
This is why we have plotted them only for W < 0.4 GeV. One has to stay below the 
two constituent quark threshold where the lack of confinement in the NJL model starts 
dominating the result. Near the threshold this effect is absent and allows for extraction 
of the polarizabilities at all tree orders. 

We show that to describe the cross sections well it is enough to consider only the first 
terms in a chiral expansion. Higher orders are important for polarizabilities. Therefore 
one can arrive at the general conclusion that the cross sections are very little sensitive 
to the polarizability. The same conclusion has been drawn in in the framework of a 
dispersive approach. Among new interesting results are also calculations of the sum of the 
electric and magnetic polarizabilities presented here both for neutral and charged pions. 
At the one-loop meson order they are equal to zero, that is why it is very important to 
advance to higher ~ p^ orders. This is possible with the NJL model, as we have shown 
here. 

Appendix. The main loop integrals 

List of the main quark-loop integrals 




(76) 




(77) 



Q4(gi, 52,93) 



h{qi,q2,q3) 




(79) 



(78) 



(80) 
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Here the notation 

Aip) 



p2 _ ^2 



has been used. A sharp euchdean cutoff A is introduced for the scalar integrals Jj, Q4, 
and all the cutoff dependence of the amplitude T^u resides in these integrals. 

We also use the standard notations for the loop functions G{s), G\{s), Ja(s) They 
are analytic in the complex s-plane, cut along the positive real axis for Re(s) > 4m^ and 
Re(s) > 4(A^ + m^) respectively. At small s, 

IQn^^^Kmy (n + l)(2n + l)! ^ ' 

le^JjU^ + m^j {n + l){2n + l)\ ^^^^ 

""^'''W^^^KA^ + my n{2n + l)\- ^^^^ 
Besides, we shall need the function 

Q3(S)=/3(P1,-P2) + — Q3(S), (85) 

s 



where 

4A2 



4(A2 + m2) - s 

References 



Ja(s) 



327r2(A2 + m2) 



+ 2 



2 

m 



A^ _|_ ffi^ 



■ (86) 



[1] J. Bijnens and F. Cornet, Nucl. Phys. B296 (1988) 557. 
[2] S. Weinberg, Physica A 96 (1979) 327. 

[3] J. Gasser and H. Leutwyler, Ann. Phys. (N.Y.) 158 (1984) 142. 

[4] J. Gasser and H. Leutwyler, Nucl. Phys. B250 (1985) 465. 

[5] Mark II Collab., J. Boyer et al., Phys. Rev. D42 (1990) 1350. 

[6] Crystal Ball Collab., H. Marsiske et al, Phys. Rev. D41 (1990) 3324. 

[7] S. Bellucci, J. Gasser and M.E. Sainio, Nucl. Phys. B423 (1994) 80. 

[8] D. Morgan and M.R. Pennington, Z. Phys. C48 (1990) 623; Phys. Lett. B272 (1991) 
134. 

[9] J.F. Donoghue, B.R. Holstein and Y.C. Lin, Phys. Rev. D37 (1988) 2423. 
J.F.Donoghue and B.R. Holstein, Phys. Rev. D48 (1993) 137. 

[10] 



A.E Kaloshin, V.V. Serebryakov, Irkutsk preprint ISU-IAP.Th 93-03, ( hep- 
ph/9306224|) . 



18 



J. Bijnens, S. Dawson and G. Valencia, Phys. Rev. D44 (1991) 3555. 

Y. Nambu and J. Jona-Lasinio, Phys. Rev. 122 (1961) 345; 124 (1961) 246. 

V. Bernard, A.A. Osipov and U.-G. MeiBner, Phys. Lett. B285 (1992) 119. 

V. Bernard, R.L. Jaffe and U.-G. MeiBner, NucL Phys. B308 (1988) 753. 

S. Khmt, M. Lutz, U. Vogl and W. Weise, NucL Phys. A516 (1990) 429; U. Vogl, M. 
Lutz and W. Weise, NucL Phys. A516 (1990) 469. 

An extensive hst of references is given in S. Klevansky, Rev. Mod. Phys. 64 (1992) 
649. 

J. Bijnens, C. Bruno and E.de Rafael, NucL Phys. B390 (1993) 501. 

J. Bijnens, preprint NORDITA 95/10 N,P ( |hep-ph/9502335|) ; J. Bijnens and J. 
Prades, Z.Phys. C64 (1994) 475. 

V. Bernard, A.H. Blin, B. Hiller, Yu.P. Ivanov, A.A. Osipov and U.-G. MeiBner, 
Preprint CRN 95-25, TK 95 18 ( |hep-ph/9506309|) to appear in Annals of Physics 
(N.Y.) (1996). 

S. Bellucci and C. Bruno, NucL Phys. B452 (1995) 626. 



A. Bel'kov, A. Lanyov and S. Scherer, Preprint JINR-E2-95-266, |hep-ph/9506i0e 
(1995). 



J. Bijnens, A. Fayyazuddin and J. Prades, preprint NORDITA - 95/84 N,P (Piep- 
ph/9512374|) . 



W. Pauh and F. Villars, Rev. Mod. Phys. 21 (1949) 434. 

V. Bernard and D. Vautherin, Phys. Rev. D40 (1989) 1615. 

C. Schiiren, E.R. Arriola and K. Goeke, NucL Phys. A547 (1992) 612. 

A. Blotz and K. Goeke, Int. J. Mod. Phys. A9 (1994) 2067; R. Lemmer and R. Tegen, 
"Meson cloud corrections to the pion electromagnetic form factor in the Nambu Jona- 
Lasinio modef, Jiilich preprint 1995. 

G. Passarino and M. Veltman, NucL Phys. B160 (1979) 151. 

V. A. Petrun'kin, Particles and Nuclei, Vol.12, Part 3 (1981) 692. 

A.E. Kaloshin, V.M. Persikov and V.V. Serebryakov, Irkutsk preprint ISU-IAP.Th 
95-01, ( |hep-ph/9504"26TD . 



B. Bajc, A.H. Bhn, B. Hiller, M.C. Nemes and M. Rosina, Z. Phys. A350 (1994) 229. 
Yu.M. Antipov et al., Z. Phys. C26 (1985) 495. 



19 



Figure captions 

Fig.l - Representative Feynman diagrams involved in the calculation of the amplitudes 
for 77 — > TTTT, subject to permutations, a) Two photons (wiggled lines) and two pions 
(dashed lines) are directly coupled to the quarks (full lines), b) Quarks rescatter with the 
quantum numbers of a scalar-isoscalar, a, in the s-channel. c) When a pair of charged 
pions is created, quarks can also rescatter in t- and u- channels to form an intermediate 
state with the quantum numbers of a pion. 

Fig. 2 - The NJL 77 tt^tt^ cross section, with averaged photon polarizations, in 
comparison with the Crystal Ball data [^. The dotted line results from the meson one- 
loop calculations (the order). The dash-dotted line includes up to tree level NJL 
corrections to the leading result. The dashed line takes into account up to p^ tree order 
NJL corrections and the full line corresponds to the total tree order corrections of NJL 
to the pion loop. 

Fig.3 - The NJL cross section for the 77 — tt+tt^ mode, with averaged photon polariza- 
tions. The dotted line corresponds to the Born approximation (p^ order), the dash-and-dot 
line up to p^, the dashed line up to p^ order. The full line denotes the complete tree order 
NJL resuh. 



20 




0.28 0.3 032 0.34 0.36 0.38 0.4 



W[GeV] 



